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We study the transport property of diffusion in a finite translationally invariant quantum sub- 
system described by a tight-binding Hamiltonian with a single energy band and interacting with 
its environment by a coupling in terms of correlation functions which are delta-correlated in space 
and time. For weak coupling, the time evolution of the subsystem density matrix is ruled by a 
quantum master equation of Lindblad type. Thanks to the invariance under spatial translations, 
we can apply the Bloch theorem to the subsystem density matrix and exactly diagonalize the time 
evolution superoperator to obtain the complete spectrum of its eigenvalues, which fully describe the 
relaxation to equilibrium. Above a critical coupling which is inversely proportional to the size of the 
subsystem, the spectrum at given wavenumber contains an isolated eigenvalue describing diffusion. 
The other eigenvalues rule the decay of the populations and quantum coherences with decay rates 
which are proportional to the intensity of the environmental noise. On the other hand, an analytical 
expression is obtained for the dispersion relation of diffusion. The diffusion coefficient is propor- 
tional to the square of the width of the energy band and inversely proportional to the intensity of 
the environmental noise because diffusion results from the perturbation of quantum tunneling by 
the environmental fluctuations in this model. Diffusion disappears below the critical coupling. 



I. INTRODUCTION 

The diffusion of particles in a condensed phase is a fun- 
damental transport process which is ubiquitous in natu- 
ral phenomena. Since the pioneering work by Einstein 
in 1905, it is known that diffusion is related to conduc- 
tion or mobility. Although diffusion has been extensively 
studied in classical systems, much less is known in quan- 
tum systems where the quantum effects can deeply affect 
the transport properties at low temperature. Yet, the 
theoretical understanding of quantum diffusion remains 
sparse and difficult because diffusion is an incoherent pro- 
cess very remote from the very coherent basic quantum 
dynamics. 

The purpose of the present paper is to study a simple 
model of quantum diffusion in a system which possesses 
the property of being invariant under spatial transla- 
tions. This property is a fundamental feature of systems 
sustaining transport processes such as diffusion or con- 
duction. The translational invariance is at the basis of 
the early studies of electronic conduction based on the 
Boltzmann-Lorentz equation 0, and is also important 
in the polaron model 0. Motivated by the problem of 
dissipation in quantum macroscopic phenomena, models 
of quantum Brownian motion have been proposed and 
studied since the eighties Q. In such models, one degree 
of freedom has a diffusive-like motion in an external po- 
tential which breaks the translational invariance. These 
models were later extended to the study of transport in 

spatially periodic potentials HEHSHHIISIlilll- 

In these studies, attention was focused essentially on dif- 
fusion. In the present study, we intend to study the 
complete set of the relaxation modes in the system by 
obtaining all the eigenvalues of the time evolution super- 
operator. The complete diagonalization of the time evo- 
lution superoperator is accomplished by fully exploiting 



the translational invariance with Bloch's theorem applied 
to the density matrices. This allows us to introduce in a 
rigorous way a wavenumber and to obtain the eigenval- 
ues which gives the Liouvillian resonances as functions of 
the wavenumber. These resonances provides us with the 
characteristic times of the relaxation toward the thermo- 
dynamic equilibrium. The hydrodynamic mode of dif- 
fusion is one among the eigenstates associated with the 
Liouvillian resonances. This mode controls the long-time 
dynamics of the system. The other modes control the 
shorter time scales and are associated with decoherence. 
In the quantum model here studied, analytical expres- 
sions can be obtained in this way for the dispersion rela- 
tion of diffusion and the other resonances. 

We notice that translationally invariant models with- 
out coupling to the environment have rec ently been 
studied in such models as quantum graphs |lj| quan- 
tum multibaker maps |14|. as well as quantum periodic 
Lorentz gases [l5|. These models describe the motion 
of a quantum particle in a spatially periodic potential. 
The energy spectrum of the particle is composed of en- 
ergy bands so that the motion is ballistic on long-time 
scales. If the energy of the particle is distributed over 
many energy bands, the motion becomes semiclassical 
and diffusion may manifest itself as a transient behavior 
before the long-time ballistic motion. This has been re- 
markably demonstrated for the quantum multibaker map 
in Ref. In order for diffusion to persist up to ar- 

bitrarily long times, the quantum system must contain 
infinitely many degrees of freedom, for instance, by the 
coupling with some environment in the subsystem-plus- 
reservoir approach we adopt in the present paper. 

The plan of the paper is the following. Our transla- 
tionally invariant model is defined in Sec. [n] by weakly 
coupling a one-dimensional tight-binding Hamiltonian to 
a delta-correlated environment. The dynamics of this 
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system is ruled by a Redfield quantum master equation 
El El III 13 HI The long-time evolution of the model 
can be studied in terms of the eigenvalues and the asso- 
ciated eigenstates of the Redfield superoperator, as ex- 
plained in Sec. IIIII The eigenvalue problem is exactly 
solved for a finite chain in Sec. IIVI The advantage of first 
taking a finite chain is that the total number of eigenval- 
ues is known and we do not miss the part of the spectrum 
controlling decoherence. The Liouvillian spectrum of the 
infinite chain is thereafter obtained in Sec. [3 by taking 
the limit N — *■ oo. Conclusions are drawn in Sec. IVfl 



II. DEFINING THE SYSTEM 



composed of the one-dimensional chain and its environ- 
ment is given by 



-fftot — H s + Hi, + A 2^ Slfy 



(5) 



where H\, is the environment Hamiltonian, Si the sub- 
system coupling operators, Bi the environment coupling 
operators, and A the coupling parameter which measures 
the intensity of the interaction between the subsystem 
and its environment. 

The dynamics of the total system is described by the 
von Neumann equation 



A. Subsystem 



dptotjt) 
dt 



AotAot(i) = -i[-fftat,/3tot(*)]) (6) 



We consider a one-dimensional subsystem described by 
the following Hamiltonian 
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represented in the site basis \l), where I takes the values 
I = 0, 1, . . . , N - 1. N is the length of the chain. We 
have here chosen periodic (Born-van Karman) boundary 
conditions. This so-called tight-binding or Hiickel Hamil- 
tonian describes a process of quantum tunneling from site 
to site and is invariant under spatial translations. 

The stationary Schrodingcr equation of the tight- 
binding Hamiltonian is given by 



H s \k) = e k \k), 



where the eigenvalues are 



e k = Eo - 2 A cos k 
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N 



and the eigenvectors 
(l\k) 



N 



(2) 



(3) 



(4) 



with k = 0,1, N — 1. Accordingly, the Hamiltonian 
l[TJl has an energy spectrum with a single energy band of 
width AA and the motion of the particle would be purely 
ballistic without coupling to a fluctuating environment. 



B. Coupling to the environment 

Now, we suppose that the subsystem is coupled to a 
large environment. The Hamiltonian of the total system 



where £tot is the Liouvillian superoperator of the total 
system. We adopt the convention that H — 1. The re- 
duced dynamics for the density matrix p(t) = TrbPtot(i) 
of the subsystem is known to obey a Redfield quantum 
master equation for weak coupling to the environment 
[H LLJ HI IH Ulj • This equation can be systematically 
derived from the complete von Neumann equation for the 
total system (0 by second-order perturbation theory so 
that 

^ = Tr b A ot p t otW h e d(t)p(t) + 0(A 3 ). (7) 
dt 

On time scales longer than the correlation time of the 
environment, the Redfield quantum master equation is 
Markovian and reads 



dp 



i[H B ,p]+\ 2 J2{Tip8l 



S\pf} - Stip - pf}s\) + 0(A 3 ), (8) 



where £r C c1 is the so-called Redfield superoperator and 
where 

POO 

T} = ]T/ drC w (T) e- i6 - T Si, e l ^ T . (9) 



The correlation function of the environment which con- 
tains all the necessary information to describe the cou- 
pling of the subsystem to its environment is given by 



C w {t) = Tr^e^B^-^Bi, 



(10) 



where is the canonical equilibrium state of the envi- 
ronment. 

The interaction of the subsystem with its environment 
is expressed in terms of the subsystem coupling operators 
which are projection operators on the site basis 



(11) 
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They take the unit value if the particle is located on 
the site I" and zero otherwise. These operators have the 
properties: 



for n = 2, 3, ... and 



sr = Si , 



£4 



(12) 



(13) 



which is the result of a coupling with delta correlation 
functions . It should also be mentioned that Eq. i|16|) 
can be directly derived from the complete von Neumann 
equation for the total system © in the singular coupling 
limit [HE!. 



III. DIAGONALIZING THE REDFIELD 
SUPEROPERATOR 



We now need to specify the environment operators by 
the choice of the environment correlation functions. We 
make two assumptions: 

Assumption 1: The environment dynamics has a very 
fast evolution compared with the subsystem dynamics. 
Therefore, the environment correlation functions decay 
in time so fast with respect to the characteristic time 
scales of the subsystem evolution, that they can be as- 
sumed to be Dirac delta distributions in time. 
Assumption 2: The environment has very short range 
spatial correlations, much shorter than the distance be- 
tween two adjacent sites of the subsystem. Accordingly, 
the environment correlation functions decay in space so 
fast that they can be assumed to be given by Kronecker 
delta in space. 

These two assumptions mean that each of the environ- 
mental fluctuations at the different subsystem sites are 
statistically independent. The environment correlation 
functions are thus given by 

C iv (t)=2Q5{t)S iv , (14) 

where Q is a real number. The operators © in the Red- 
field equation JSJ) therefore become 



Ti = QSi 



(15) 



It has been shown in Ref. 21] that the form taken by 
the operator (| 1 51) can be physically justified when r t h <C 
ib <Ct s , where Tth = ^/k^T is the thermal time, it, the 
correlation time of the environment or bath, and r s the 
subsystem characteristic time. 

Because of the fast decay of the temporal and the spa- 
tial correlations 114fl and the properties of the subsystem 
coupling operators (|llf) . the Redfield equation we have 
to solve takes the form 



dp 



= C-KedP 

= -i[H,J] 



(16) 



-\ 2 QY,(2SipS t -S, 



psf) 



0(A 3 



It can easily be verified by projecting this equation onto 
the site basis that it is translationally invariant (shifting 
all the site indices appearing in the projected equation 
by a constant does not modify the equation). Further- 
more, this equation preserves the complete positivity of 
the density matrix because it has the Lindblad form [2^ 



The eigenvalues s„ and associated eigenstates p v of the 
Redfield superoperator are defined by 



C 



Rod P 



P 



(17) 



where v is a set of parameters labeling the eigen- 
states. Because the Redfield superoperator is not anti- 
Hermitian, its eigenvalues can be complex numbers with 
a nonzero real part. The eigenvalue problem of the Red- 
field superoperator is important because the time evolu- 
tion of the quantum master equation can then be decom- 
posed onto the basis of the eigenstates as 



P(t) 



— g^Rec 



N 

t p(0)=^c,(0)e^. (18) 



The dynamics is therefore given by a linear superposition 
of exponential or oscillatory exponential functions. Since 
the reduced density matrix of the subsystem has TV 2 ele- 
ments, there is a total of N 2 eigenvalues and associated 
eigenstates. 



A. Bloch theorem 

Since the system is invariant under spatial translations, 
we can apply the Bloch theorem to the eigenstates of 
the Redfield superoperator. Thanks to this theorem, the 
state space of the superoperator can be decomposed into 
independent superoperators acting onto decoupled sec- 
tors associated with a given Bloch number, also called 
wavenumber [lj. 

Wc define the superoperator T a of the spatial transla- 
tion by a sites along the system (a is an integer) as 



(TaP V )w - P"+a,l>+a , 



(19) 



where we use the notation (Z|p|/') = pw ■ This superop- 
erator has the group property 



"~^a~^a' r ^~a f '^~a -^a-t-a' 



(20) 



Because of the translational symmetry of the system, the 
translation superoperators commute with the Redfield 
superoperator 



Red 



= . 



(21) 



4 



Therefore, the Redfield superoperator as well as the 
translation superoperators have a basis of common eigen- 
states. If r(a) denotes the eigenvalues of the translation 
superoperator, we have that 

t p" = r(o) p u , (22) 

where, because of the unitarity of T ai 

|r(a)| 2 = l. (23) 

Equation 12(J|) implies 

r(a + a ) = r(a) r(a') . (24) 

Because of Eqs. (j23> and lf2*i|l. we find that 

r(a) = e tqa , (25) 

where q is the Bloch number or wavenumber, whereupon 
we get 



Pl+a,l'+a — c Pw 

A useful consequence is that 

Pw = z iql ph-i ■ 



(26) 



(27) 



In order to determine the allowed values of the Bloch 
number, we write by using Eq. that 



fe.fe' 



i(lk-l'k')'- 



and 



Pi+\,i'+\ 



fe,fc' 



i(lk-i'k')^- e i(k-k') 



Because of Eq. ()26[) . we also have 



(28) 



■(29) 



(30) 



Multiplying both sides of Eqs. and (39 by 

(l'\k"')(k"\l), taking the sum J^i v °f it, and identifying 
them, we obtain 



l9 (fc|/5Hfc') 



a i{k-k')i 



\p»\k<). 



(31) 



We can now notice that if q ^ (k-k')jj-, then {k\p v \k') = 
0. Finally, using the periodicity 



Pl+N,V — PlV i 

Pli'+n — Pw > 



(32) 
(33) 



and Eq. I|26|) . we find that the Bloch number takes the 
values 

q = jjf, where j = 0, 1, . . . , iV - 1 . (34) 

Consequently, the Redfield superoperator can be block- 
diagonalized into N independent blocks, which each con- 
tains N eigenvalues as we shall see in the following. 



B. Simplifying the problem 

The eigenvalue problem of the Redfield superopera- 
tor can be formulated in each sector labeled by a given 
wavenumber q. For this purpose, Eq. I|17|) with the ex- 
plicit expression i|16|) of the Redfield superoperator is pro- 
jected onto the site basis to get 

SuPw = ~iA {-pi-. x ,v - Pi+i,v + Ph'-i + Pu'+i) 

+2\ 2 Q{8 w -l) Pw ■ (35) 

Using Eq. I|27[) and replacing V — I by I, we have 

( S „ + 2QA 2 )^ = 24 (sin |) (e"^ +1 _ e^p^i-i) 
+2Q\ 2 p» 5 0l . (36) 
Making the change of variable 

Poi = r'e^'/i , (37) 
we obtain the simpler eigenvalue equation 

Pvfi = Soifi-iP(fi-i+fi+i) , (38) 

where 



p v 



2QA 2 



1 



and 



, A . q 
= QA 2 Sm 2 ' 



(39) 



(40) 



IV. FINITE CHAIN 

A. The eigenvalue problem 

The expression (|38ll can be written in matrix form 
without the index v to simplify the notation, 



pf = Wf, 



(41) 



where p denotes the eigenvalue, / = (/o, . . . , /jv-i) the 
eigenvector of size N, and W the following N x N matrix 



/ 1 -i(3 

-i(3 -i/3 

-0 -0 



-0i~ N e iN i \ 



-0 -0 
-0 
-0 



-0 



(42) 



We look for eigenstates of the form 
fl = Ac m + Be~ m 
Solving Eq. gJJ with gSJl gives 



(43) 
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• for < I < JV - 1: 

/i = — 2i (3 cos 9 

• for Z = 0: 



(44) 



(1 - n)(A + B) - i(3{Ae lB + Be- 10 ) 
-t(3t- N e m i(Ae^ N -^ + Be~^ N -^) = , (45) 

• for Z = AT — 1: 

Ae * eN + Be-* eN - i N e~ iN % (A + B) = . (46) 

Solving the homogeneous linear system of Eqs. (|45|l and 
and replacing /j by one gets the characteristic 
equation 

2i/3 sin [R{q) - cos 6>JV] = sin 9N, (47) 

with 

*(«) = \ ( 
Using Eq. ijSljl. we find 



■iV ~iNZ , ~A r iJV- 

« e 2 + z e : 



(48) 



for JVodd: i%) = , (49) 
for N = 41: R(qj) = , (50) 

for N = 4I + 2: R( 9j ) = , (51) 

with / integer. From now on, we shall speak of even 
(respectively odd) q, if q corresponds to an even (respec- 
tively odd) integer j in Eq. 1|34[> . Therefore, the charac- 
teristic equation (|4*7I) becomes 



for TV odd: 



B. The diffusive eigenvalue fi^ 

We first look for an eigenvalue /i which is real and 
should correspond to a monotonic exponential decay. 
With this goal, we suppose that the angle 6 is complex 

= Z + iri, (57) 

so that the eigenvalue l|44|l becomes 

fi = —2if3cos9 — —2(3 sin £ sinh?7 — 2«/3cos£ coshr; .(58) 

This eigenvalue is real under the condition that cos £ = 
which is satisfied for 



in which case 



(i = —2 (3 sinh?; 



(59) 



(60) 



We notice that the condition cos £ = is also satisfied for 
£ = § but it can be shown that this other case leads to the 
same eigenvalue as (|59|) . If we introduce the conditions 
(|57jl and |J5S) in Eqs. (52J|, lf53)l. and lj53)l. we get 

• for JV odd: 

2/3 cosh 77 = coth/V^ ; (61) 

• for N — 41 and g even or for N = 41 + 2 and q odd: 



2/3 cosh 7/ = coth 



iV?7 



(62) 



• for N = 41 and q odd or for JV = 41 + 2 and q even: 



2i/3sin0 = -tan6W 



(52) 



2(3 cosh 77 = tanh 



JV77 



(63) 



for JV = 41 and g even or for JV = 41 + 2 and g odd: 
either 



or 



6W 

2idshi9 — cotan 

2 



;6»JV = 1 



(53) 



(54) 



• for JV = 41 and q odd or for JV = 4/ + 2 and q even: 
either 



6»JV 

2i(3sm6 = —tan 

2 



or 



cos ON = -1 . 
We solve Eqs. (|521l - (|5o]) as follows. 



(55) 



(56) 



In the limit JV — > 00, the right-hand side of these equa- 
tions tends to unity if a nonvanishing solution 77 7^ ex- 
ists. In this limit, this solution is thus given by 



arccosh — , 

2(3 



which exists only if (3 < |. Because 



sinh 770 = 



-1, 



the corresponding eigenvalue should be 
Mo = Vl - (W ■ 



(64) 



(65) 



(66) 



However, for a finite chain with JV < 00, we expect 
a correction <5?7 to the solution 77 = 7/0 + Srj. Replacing 
this correction in Eqs. I|61|) . Q62p. and l|63|) . we obtain by 
Taylor expansion that 
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• for N odd: 



Srj 



^— 2iVarccosh ^ 



(67) 



• for N — 41 and q even or for TV = 47 + 2 and q odd: 

„— iVarccosh 



<5r?~2 



(68) 



• for N = 41 and g odd or for N = 47 + 2 and g even: 

— TVarccosh ^ 

5r) ~ -2 -^=_=- ; (69) 

up to corrections of 0(<5r; 2 ). Using the expression l)6()[l. 
we finally obtain the eigenvalue 

/x (1) = 2/3 sinh 770 + 2/3 cosh ry <5?7 + 0(^ 2 ) 

= v/1 - (2/3) 2 + <fy + 0(<5r/ 2 ) , (70) 

where (5r/ is respectively given by Eqs. (|67|l . (|68|l . and 
(I69|l . Accordingly, the correction 5r] to the eigenvalue 
decreases exponentially fast with the size N of the chain. 

Using Eqs. I|39[) and 1)40 [I. we finally obtain the eigen- 
value 

,.0) = _ t .. 

/ / 1 /I „\ 2 

- 2QA 2 + 0(«5t7) . 

(71) 

The eigenvalue = corresponding to a vanishing 
wavenumber q = is always in the spectrum of the Red- 
field superoperator. The associated eigenstate describes 
the stationary equilibrium state. At low wavenumbers 
q, (3 — > 0, we recover the dispersion relation of diffusion 




«W = -Dq 2 + 0(q 4 ) , 
with the diffusion coefficient 

D A2 



(72) 



QA 2 



which justifies calling fi^ or the diffusive eigenvalue. 

We notice that the diffusive eigenvalue no longer exists 
beyond the critical value f3 c = h. Since the matrix H42|) 
has a total of N eigenvalues, we expect further nondiffu- 
sive eigenvalues in a number of N — 1 for (3 < \ and N 
for (3 > |j as confirmed in the following subsections. 



C. The eigenvalues /v' 

Beside the diffusive eigenvalue, we expect eigenvalues 
corresponding to the chain-like structure of the matrix 



(|42|l . To obtain these eigenvalues, we pose tan% = 
2?/3sin0 so that Eqs. lj52*|) . and can be written 
respectively 



sin(0A + x) = , 

/6W A 
cos I — + x I =0 

/0iV A 
sm I — + X I = . 



(74) 
(75) 

(76) 



7T 


9N 


+ 2" 






ON 


n7r — 









Therefore, we obtain 

i(nir-9N) = arctanh(-2/3sin0) , (77) 

/ j //;;- [ - - — ] = arctanh(-2/3sin0) , (78) 

= arctanh(-2/3sin0) . (79) 

We now expand arctanh(— 2/3 sin 0) around (3 = 0: 

arctanh(-2/3sin0) =° -2/3 sin - ^/3 3 sin 3 + 0((3 5 ) . 

(80) 

If (3 = 0, the solutions of Eqs. lff7|l. lf?8}. and (O are 

respectively given by 



7,| : — — . \\ IK'i <' // = 1, Z, . . . , JV — 1, (81) 



— , where n=l,2 N — 1, 

(2n+l)7T . A" 



v , where n = 0, 1, 1, (82) 

2mr N 
9 " = ~N~' Where n = 1 ' 2 '---'"2~ _1 ' (83) 

Notice that n = is rejected in Eqs. (HU and JHSJ- It 
is due to the fact that = does not correspond to an 
eigenvector because it can be seen that /; = A + B ^ 
in Eq. I|43(l cannot be an eigenvector of Eq. 1421) . 

Using the expansion (fHOJ in Eqs. fTTjl. f7B|l and 1(79} 
with = O + 89, we find 

/3^o i2^ n i8f3 3 3 



sm f n sin t/n 

M 3M 



(84) 



where M = N for Eq. J77J and M = f for Eqs. (J7HJ) and 
(73) (|79|l . The eigenvalue 144(1 is now given by the expansion 



H = -2i/3cos(0 o + 69) 

S9 =° -20 cos O + 20 sin O 50 

+i/3cos0 o S9 2 + 0{59 3 ). (85) 
Using Eq. (g3J in (JHSJl gives 



^(2) /3^o _2i/3cos0 o + [ i 



M 2 

3 1 
3M 



sin 2 O + ^ sin 4 O f O ( - 



(86) 



Consequently, we have 



7 



• for JV odd, using Eq. {ggj with (fHT|) : 



E. The eigenvalues /i 



(4) 



M (2) _ 2 - /3cos 



7J7T 4/3 2 . 2 7T,7T 16/3 4 . 4 7T.7T 



JV JV JV 3JV JV 
where n = 1, 2, . . . , JV - 1 ; (87) 

• for N = 41 and (? even or for JV = 4/ + 2 and g 
odd, using Eq. (JSgJl with (JH2J): 

( 2 > -2z/3cos (2n+1) " 



JV 



8/3 2 . 2 (2n+l)7r 32/3 4 . 4 (2n + lW 
AT Sm -^V— + W Sm ^V— ' 



JV 

where n = 0, 1, . . . , — — 1 



(88) 



• for JV = 4/ and q odd or for JV = 41 + 2 and q even, 
using Eq. JggJ with JggJ: 

(2) /3-+o „.„ 2n7r 8/3 2 . 2 2mr 32/3 4 . 4 2mr 

u K > = -2i/3cos h -— sm 1 — sin 4 , 

M JV JV JV 3JV JV' 



JV 

where n = 1, 2, . . . , — — 1 . 



(89) 



An important observation is that, for /3 < i, the ex- 
pansion H80fl which implies 2/3 sin < 1 is satisfied every- 
where, i.e., for all the values of 9 and therefore for all the 
eigenvalues. However, for j3 > i, the Taylor expansion 
around /3 = in Eq. JggJ is only valid if 2/3 sin 9 < 1 . 
Therefore, a transition zone exists around sinff ~ 1/(2/3). 
According to Eq. (|4*4fl . this transition corresponds to the 
critical value of the eigenvalue given by 



Me = ±V(2/?) 2 - 1 



(94) 



Therefore, for (3 > |, the expansion l|8U|) around /3 = is 
only valid if \fj,\ > |/i c |. For < |// c |, we should instead 
consider the asymptotic expansion of arctanh(— 2/3 sin 0) 
around (3 = oo: 

arctanh(— 2/3 sin 9) ^=°° /_ _ 



2 2/3 sin (9 24/3 3 sin 3 ' 

1 °(w 



160/3 5 sin 5 i 



(95) 



(3) 



D. The eigenvalues /i 



The solutions of Eq. I|54|l are simply given by 

2mr JV , x 

0=—, where n = 1, 2, 1 . (90) 

We reject n = and n = y because the corresponding 
eigenvector does not exist in these cases. Similarly, the 
solutions of Eq. (|56|l are given by 



which leads to another family of eigenvalues existing for 
/?>§. 

If j3 -> oo, the solutions of Eqs. if??}. (JZgJ and (JZHl are 

respectively given by 



(n + |)tt 



JV 



where n = 0, 1, . . . , JV - 1 , (96) 



2nir N 
7 = -j^-, where n = 1, 2, - 1 , (97) 



JV 

(2n+ 1)tt 

JV : 



JV 

where n = 0, 1, . . . , — - 1. (98) 



„ (2n + l)7r JV , x 

g= N , where n = 0, 1, 1 . (91) 

Consequently, we have the further eigenvalues: 



• for JV = 41 and g even or for JV = 4/ + 2 and q 
odd, using Eq. (gtj with i|5U|): 



/' 



2n7r 
1\T 



(3) = -2i/?cos [ ) , where n = 1, 2, - 1 



JV 

T 



(92) 



Because of the condition \fi\ < \/i c \ with the critical values 
l)94[l . we should only consider the angles in the interval 
$o,c < 9q < n — 6*0,0 with 



70, C 



arcsm ■ 



1 



(99) 



so that the integer n in Eqs. JggJ, (|9"7|) . and lj9"51) is 

restricted to take the intermediate values n m ; n < n < 
n max which do not reach the extreme values. 

Using the expansion @2J in Eqs. ijTTjl. (JJgJ, and JZHJ 
with = 6*o + (50, we find 



• for N — 41 and q odd or for JV = 4/ + 2 and q even, 
using Eq. JSJ with (fSTfl : 



.(3, = _^ C03 ((5l±i>), 



JV 

where n = 0, 1, . . . , — — 1 . (93) 



2/3Msin6» 24/3 3 M sin 3 6> 160/3 5 Msin 5 6 
i \ „ / 1 



-0 



^ 7 Jlf 



/3 2 m 2 y ' 



(100) 



where M = N for Eq. G7J and M = § for Eqs. O 
and CgJ. 
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Using Eq. 1 100(1 in (JHSJ gives 



Inserting in Eq. (|85|l . we obtain 



(4) 0- 



-2if3cos9 + — ( 1 



O 



12/3 2 sin 2 6» 



4 sin 4 9 



(3 6 M 



O 



(3M 2 , 
(101) 



Consequently, we have 

• for N odd, using Eq. (|TUT)l with 



/./, 



( 4 ) /3^oo 



2z/?cos 



(n+ i)7T 




1 



12/? 2 sin 2 



2V 



where n = 0, 1, . . . 



,iV-l; 



(102) 



• for N ~ 41 and g even or for AT = 4/ + 2 and g 
odd, using Eq. ((TUTJ) with 

(4) /3^oo 2717T 

/ w = — 2ipcos 




N y ' 12/3 2 sin 2 2j™L 
N 

where n = 1, 2, . . . , — — 1; 



4 sin 4 2^ 



(103) 



• for N = 41 and q odd or for TV = 41 + 2 and q even, 
using Eq. l(TUT|) with 

«W ^ -2i/3cos^±^ 



iV 



A/ 



l 2/3 2 sin 2 



80/34 sin 4 (2"+l)£ 



AT 



where n = 0,1, '— - 1 ; (104) 

with the aforementioned restriction on the values of the 
integer n. 



F. The eigenvalues 

For N = 41 and q even or N = 41 + 2 and q odd, 
two special eigenvalues exist in the limit /3 — > oo around 
6>o = and #o = 7r. They can be obtained by taking 
= 8 + 58 and directly solving Eq. to get in both 

cases 



(105) 



G. Description of the spectrum 

Rewriting the eigenvalues l|4~4*)l of the Redfield super- 
operator with their explicit dependence in terms of Eq. 
||tnj|, we get 



2i(3 cos = -2i—, r (sinf) cos (9 , (107) 
QX V 2/ 



where v can take A^ 2 different values because q and 9 
take AT values each. Remembering that according to Eq. 
(EH 



2QA 2 (av - 1) 



(108) 



we conclude that we have found in this section all the 
eigenvalues of the Redfield superoperator. We list them 
in Tables [fl and ^ according to the parameter regime in 
which they hold. 

We now discuss the main features of the spectrum 
when the different physical parameters are varied. This 
discussion is based on our analytical results for the eigen- 
values and on the comparison between these results and 
the eigenvalues obtained by numerical diagonalization of 
the Redfield superoperator. Since the eigenvalues /i q g are 
related to the Redfield superoperator eigenvalues s q e by 
Eq. (|108fl . we notice that all the eigenvalues of the com- 
plete spectrum always satisfy < Re fi q g < 1 or, equiv- 
alently, — 2QA 2 < Re s q g < 0. The imaginary part of s q g 
is simply proportional by a factor 2QX 2 to the imaginary 
part of fi q g . 

We start by studying the N eigenvalues \i q e obtained 
by fixing the wavenumber q (even or odd) and varying 8. 
For given physical parameters (A, A, Q, N), fixing q is 
equivalent to fixing (3. 

For /3 < \ , the analytical expressions of the eigenvalues 
which concern us are summarized in Table [I] Two fami- 
lies of eigenvalues (/i^ and enter in the discussion 
for N odd, and three families (p^ 1 ', p 2 ', and ^ 3 ' ) ), for 
N even. The numerical eigenvalues are plotted in Figs. 
ma),[2fa), and|3Ja) and are in very good agreement with 
the analytical results. The sole diffusive eigenvalue /z^ 
has a real part and no imaginary part. The AT — 1 other 
eigenvalues, either belongs to the ijS 2 ' family for N odd 
or to the /i' 2 ' and fj,^ families for N even. The eigenval- 
ues and /j,^ have an imaginary part which extends 
from —2/3 to 2(3 and they generate oscillations in the dy- 
namics. The real part of the eigenvalues is small and 
tends to zero in the large N limit. The real part of the 
M ( 3 ) ei genvalues is always zero. 

For /3 > i, the diffusive eigenvalue /Z 1 ) has disap- 
peared after merging with the other eigenvalues and the 
situation is slightly more complicated. The situation for 
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a moderate value of (3 > | is depicted in Figs, ^b), 
I2{b), and Elb) while the analytical expressions of the 
eigenvalues are given in Table [H] Since fi^ no longer 
exists, we have the two families of eigenvalues fj,( 2 > and 
^( 4 ) if TV i s odd, and the three families ^ 2 \ /i*- 3 -* and 
^(4) if i s even. Two regions of the spectrum have to 
be distinguished. The eigenvalues /i*- 2 ' exist in the region 
where > |/i c | while the eigenvalues ^ 4 ' exist in the 
region where < |/i c |. We observe that the extra family 
of eigenvalues // 4 ) has appeared because of the collision 
with the diffusive eigenvalue fjP-' . We can see in Figs. 
^b),[2Ib), and^b) that the analytical results of Table 
[njreproduce very well the eigenvalues obtained by numer- 
ical diagonalization in the two regions. Here again, the 
number of eigenvalues is equal to N for a given wavenum- 
ber q, the imaginary part of the eigenvalues extends from 
—2/3 to 2/3, and the real parts of all eigenvalues tends to 
zero in the large N limit. 

A special situation occurs when (3 > ^ is increased to 
large values. This situation is depicted in Figs, ^c), 
|2Jc), and[3{c). The situation is similar to the previous 
one but the region > |/i c | has disappeared so that 
the family of eigenvalues /i*- 2 -* corresponding to the ex- 
pansion j3 — *■ no longer exists. For N odd and for N 
even with q odd, these eigenvalues are replaced by the 
eigenvalues fi^ eigenvalues. For N even and q even, we 
find the two eigenvalues fi^ beside the family of eigen- 
values /i^ 4 -* . The agreement between the analytical and 
numerical results is very good here also. As before, the 
imaginary part of all these eigenvalues extends from — 2i/3 
to 2i(3 and their real parts tends to zero in the large N 
limit. 

A global view of the complete spectrum of the N 2 
eigenvalues of the Redfield superoperator is depicted in 
Fig. 01 by varying the wavenumber q in a third dimen- 
sion. Here, we only consider for simplicity the case where 
N is odd. The relation between the wavenumber q and 
the parameter /3 is given by Eq. I|4l)[) . The wavenum- 
ber q varies in the first Brillouin zone or, equivalently, in 
the interval < q < 2tt. We see in Fig. 0Ja) that the 
diffusive eigenvalues fi^ exists for all the values of the 
wavenumber in the case < \ which implies (3 < \. 

However, if > |, the diffusive eigenvalue dis- 
appears as expected for some values of the wavenumber 
corresponding to j3 > i. This situation is observed in 
Fig. lib). 

For very large values of > i, the diffusive branch 
of the spectrum is reduced to the sole eigenvalue at 
gW = 0, as seen in Fig. 01c). In this case, diffusion 
has disappeared from the spectrum which only contains 
eigenvalues associated with damped oscillatory behav- 
ior. The diffusive branch can be supposed to have disap- 
peared when its last nonzero eigenvalue disappears in Eq. 
(|Tip . Therefore the diffusive branch disappears when the 
value of (3 for the first nonzero eigenvalue correspond- 
ing to q = ^ is larger than the critical value /3 C = \. 
This happens when the coupling parameter exceeds the 



critical value given by 

\2A tt N>b \2A-k 

A ==V^ sm » " vw- <109) 

This disappearance of the diffusion branch can be ob- 
served in Fig. 01c). We notice that the diffusive branch 
always exists in the infinite-system limit (N — > oo) in 
which case A c can be arbitrarily small. 

We have described in this section the complete spec- 
trum of the Redfield superoperator for a finite chain. We 
now briefly indicate which are the dynamical implications 
of these results. 



H. From the spectrum to the dynamics 

The linear decomposition l|18|) of the density matrix 
shows that the modes which control the long-time dy- 
namics correspond to the eigenvalues having the small- 
est absolute value of their real part. The eigenvalues 
with larger absolute value of their real part correspond to 
faster modes which control the relaxation on shorter time 
scales. For systems of finite size N, the long-time relax- 
ation can be of two kinds. In the case where A < A c , the 
long-time relaxation is controlled by nondiffusive modes 
and consists in complicated oscillations of different peri- 
ods damped at rates Res ~ —2QX 2 . This is due to the 
fact that, in this case, all the eigenvalues of the spectrum 
have a similar real part and different imaginary parts. 
This also indicates that the modes corresponding to the 
relaxation of the coherences (modes with complex eigen- 
values) and of the populations (modes with real eigenval- 
ues) decay on similar time scales. In the other case where 
A > A c , the long-time relaxation is controlled by the diffu- 
sive mode. This relaxation is free of any oscillations and 
is controlled by the rate s ~ ~<in 2 A 2 /(Q\ 2 N 2 ). This 
diffusive relaxation exclusively concerns the populations 
of the system. The other nondiffusive modes describe the 
decoherence (as well as the relaxation of the populations 
beside the part controlled by diffusion when 2 A > QX 2 ). 
These modes decay at rates Re s ~ —2Q\ 2 , i.e., on much 
shorter time scales than the diffusive mode. The dynam- 
ics of the finite N system is studied in detail elsewhere 

m 



V. INFINITE CHAIN 

The spectrum of the infinite chain coupled to its envi- 
ronment can be obtained from the spectrum of the finite 
chain in the infinite-size limit N — ► oo. The wavenumber 
q becomes a continuous parameter varying in the first 
Brillouin zone — ir < q < tt. 

For given wavenumber, the diffusive eigenvalue (j,^ or 
S W given by Eq. (|71|l remains isolated. Consequently, we 
obtain the result that the dispersion relation of diffusion 
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is exactly given by the analytical expression 



2yQ 2 A 4 - (2Asin|) -2QX 2 = -Dq 2 + 0(q 4 ) . 

(110) 



The diffusion coefficient 



D 



A 2 

QX 2 ' 



(111) 



is proportional to the square of the parameter A of the 
tight-binding Hamiltonian and inversely proportional to 
the parameter QX 2 of the coupling to the environment. 
The transport is therefore due to the tunneling from site 
to site, which is hindered by the environmental fluctu- 
ations proportional to QX 2 . It as been shown in Ref. 
|2l| that, for an Ohmic coupling to the environment, the 
diffusion coefficient is inversely proportional to the tem- 
perature. By using the Einstein relation between the 
diffusion coefficient and the conductivity, this latter is 
therefore inversely proportional to the square of the tem- 
perature. 

For 2 A < QX 2 , the diffusive eigenvalue exists for all 
the values of the wavenumber — n < q < +n, as seen in 
Fig. Eta). 

For 2 A > QX 2 , the diffusive eigenvalue only exists for 
all the values of the wavenumber in the range 



< a < 



-'7c 



with q c — 2 arcsin 



QX 2 
2A 



(112) 



as seen in Fig. OJb). 

Beside the isolated diffusive eigenvalue, the spectrum 
at given wavenumber q contains a continuous part ob- 
tained by the accumulation of the eigenvalues ^ 2 \ fi^> , 
fj,^', and /i^ in the limit N — > oo. Indeed, in this 
limit, all these eigenvalues accumulate into a segment 
of straight line given by 



-MA 



( Sin |) 



cos(9-2QA 2 



(113) 



with < 9 < tt and -7r < q < +ir [see Eqs. (fTTTTf) 
and 1|108|) ]. This part of the spectrum is also depicted 
in Fig. El and describes the time evolution of the quan- 
tum coherences which are damped at the exponential rate 
—Re s q g — 2QX 2 with possible oscillations due to their 
nonvanishing imaginary part Im s q g. 

The eigenvalue 1110|l can be considered as a Liouvil- 
lian resonance [25l l2fil] similar to the Pollicott-Ruelle res- 
onances describing diffusion in classical systems [221 123 ■ 



VI. CONCLUSIONS 

In the present paper, we have studied an exactly solv- 
able model of simple translationally invariant subsystems 
interacting with their environment. The coupling to the 
environment is described by correlation functions which 



are delta-correlated in space and time. The reduced dy- 
namics of the subsystem is described by a Redfield quan- 
tum master equation which takes, for such environments, 
a Lindblad form. Thanks to the invariance under spatial 
translations, we can apply the Bloch theorem to the sub- 
system density matrix. In this way, we succeeded in get- 
ting analytical expressions for all the eigenvalues of the 
Redfield superoperator. These eigenvalues control the 
time evolution of the subsystem and its relaxation to the 
thermodynamic equilibrium. Two kinds of eigenvalues 
were obtained: the isolated eigenvalue (|110fl giving the 
dispersion relation of diffusion along the one-dimensional 
subsystem and the other eigenvalues (|113fl which describe 
the decay of the populations and quantum coherences. 
The process of decoherence in the subsystem is controlled 
by these latter eigenvalues l|113|) . 

The properties of the system depend on the length N of 
the one-dimensional chain, on the width AA of the energy 
band of the unperturbed tight-binding Hamiltonian and 
on the intensity Q of the environmental noise multiplied 
in the combination QA 2 with the square of the coupling 
parameter A of perturbation theory. 

We discovered that, for a finite chain, there are two 
regimes depending on the the chain length N and the 
physical parameters A and QA 2 . 

For a finite and small enough chain, there is a non- 
diffusive regime characterized by a time evolution with 
oscillations damped by decay rates proportional to QX 2 . 
The oscillations are the time evolution of the quantum 
coherences. This nondiffusive regime exists if the cou- 
pling parameter is smaller than a critical value which is 
inversely proportional to the square root of the chain size 
N: X < A c = O(N-i). 

For larger chains, we are in the diffusive regime with 
a monotonic decay on long times at a rate controlled 
by the diffusion coefficient. In this regime, the slower 
relaxation mode relaxes exponentially in time with the 
scaling t/(XN) 2 . 

In the limit of an infinite chain N — > oo and for 
non- vanishing coupling parameter QX 2 , the nondiffusive 
regime disappears and the system always diffuses. 

The diffusion coefficient is proportional to the square 
of the width 4A of the energy band and inversely propor- 
tional to the intensity QX 2 of the environmental noise. 
Accordingly, we are in the presence of a mechanism of 
diffusion in which the quantum tunneling of the parti- 
cle from site to site is perturbed by the environmental 
fluctuations. 

The eigenvalues of the Redfield superoperator obtained 
in the present paper give the Liouvillian resonances at 
the second order of perturbation theory. In this reg ard, 
the present work extends the results of Refs. |25l l2fil ] 
on the spin-boson model to systems with a translational 
symmetry in space and capable of sustaining the trans- 
port property of diffusion beside simple decay processes. 
These Liouvillian resonances are the quantum analogues 
of the Pollicott-Ruelle resonances which have been stud- 
ied elsewhere for diffusion in classical systems [271 H3 • 
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TABLE I: For wavenumbers q corresponding to j3 < |: List of the eigenvalues fi q e of the matrix II42II given by Eq. (111)71 . 

These eigenvalues are directly related to the Redfield superoperator eigenvalues by Eq. 110811 . 



TV odd 



TV = 41 and q even 

or TV = 41 + 2 and q odd 



TV = 41 and q odd 

or TV = 41 + 2 and q even 



M (1) = V l- (2/3)2 + 2 *-~ Bhl/ 



Vl-(2/?) 2 + 2^= 



sh 1/(2/3) 



(2/3)2 



= Vl~(2/3) 2 



-Narccosh l/(2.rf) 



Vi-(2/3) 2 



,(2) ^o_ 2i/ 3 cos - 

+ Mi S i n 2 2ZL 

T N SID 

3 

+0(#)+0(f^) 



for n = 1, 2, 



,/V- 1 



/ 2 > ^=°-2i/3cos^±i^ 

, 8/3 2 -2 (2n+l)7r 
T sin jy 

32/3 4 o;ri 4 (2n+l)7r 



for n 



AT 



+ 0(fr) 
0,1,..., f 



,( 2 ) Pz 



-2i/3 cos 2^ 



8j3- 



N 
^4 2n7r 



to('f) + o(^ 

for n = 1, 2, . . . . 



2 



/i (3) = -2i/3cos(2^L) 
for n = 1,2, ...,f 



= -2 l /3cos((22±l^) 
for n = 0, 1, . . . , f - 1 



TABLE II: For wavenumbers q corresponding to j3 > |: List of the eigenvalues (j, q e of the matrix 11 2t given by Eq. I|l()7^ . 

These eigenvalues are directly related to the Redfield superoperator eigenvalues by Eq. HUSH . 



TV odd 



TV = 47 and q even 

or TV = 41 + 2 and q odd 



TV = 41 and g odd 

or TV = 41 + 2 and <? even 



If | M (4) < |Mc|:° 



M (4) 3-oc_ 2 .^ cos (n+4), 



12/3 2 sin 



32/3 4 sin' 



(n + (l/2))i. 



+0(^) + 0(^) 
for n min < n < n max 



,(4) /3" 



-2i/3cos2^L 



+#(1 + 



12/32 sin 2 2^ 



J 1 

~^32/3 4 sin 4 2jpE' 

+0(^) + 0(^) 
for n min < n < n max 

M ( 5 ) ^^+1+0(1) 



,(4) 0- 



_2j/3cosi22±l)iL 



12/3 2 sin 



2 (2n+l)n 



32,34 sin 4 

+0(^) + 0(^) 
for n min < n < n max 



If | M (2) > |Mc|: 



^ =°-2if3cos^ 



sin 2 



I 4/3- 
, 16/3 4 



+0(^)+ (^ 

for n = 1, 2, 



and n — n n 



.,TV-1 



/■' 



(2) /3=r° 



- -2i/3cos (2 "+ 1)7r 

, 8/3 2 -2 (2n+l)7r 
~r AT hlU JV 
I 32/3 4 -4 (2n+l)7r 
3AT „ S111 JV, 

for n = 0, 1, . . . , n m in 
and n = n max x 



•' 2 



— - 1 



(2) /3^° o-o 

= — 2i/3cos 



2 » ~ 
2V 

+ Ml S i n 2 a»* 



N 

+ f£sin 4 2 ' ' 



+0(f) + 0(f^) 
for n = 1, 2, 



and n = n n 



, ^min 

N 



— — 1 



/ 3) = -2i/3cos(^) 



for n = 1, 2, . 



' 2 



M < 3 > =-2i/3cos(i2n+iM) 



for n = 0, 1, 
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FIG. 1: Eigenvalue spectrum for TV = 21 and given q: (a) 
/3 = 0.4, (b) /3 = 0.7, and (c) j3 = 10. ju e denotes the exact 
eigenvalues obtained by numerical diagonalization and //^ 
the eigenvalues of the different families given in Tables □ and 
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FIG. 3: Eigenvalue spectrum for N — 20 and q odd: (a) 
/3 = 0.4, (b) P = 0.7, and (c) j3 = 10. fi c denotes the exact 
eigenvalues obtained by numerical diagonalization and jit' 1 ' 
the eigenvalues of the different families given in Tables Eland 

El 




FIG. 4: Complete spectrum of the eigenvalues s q e versus the 
wavenumber q for the Redfield superoperator 1161 - 117( 1 The 
eigenvalues are given by Eq. 110811 in terms of the eigenvalues 
fi q e depicted in Figs. Ill3l at given values of the wavenumber. 
The size of the chain is here N = 21. The parameter values 
are: (a) ^ = 0.4; (b) ^ = 0.7; (c) ^ = 10. 



FIG. 5: Spectrum of the infinite chain coupled to its envi- 
ronment: (a) in the regime 2A < QX 2 for A — 0.4 and 
QX 2 = 1; (b) in the regime 2A > QX 2 for A = 0.6 and 
QX 2 = 1 where the diffusive branch is limited to the wavenum- 
bers \q\ < q c = 1.97022. 



